Revising neutrino oscillation parameter space with direct flavor-changing interactions 
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We formulate direct, neutrino flavor-changing interactions in a framework that fits smoothly 
with the parameterization of two- and three-state mixing of massive neutrino states. We show that 
even small direct interaction strengths could have important consequences for the interpretation of 
currently running and proposed oscillation experiments. The oscillation amplitude and the borders 
of the allowed regions in two- and three-flavor mixing parameter space can be sensitive to the 
presence of direct interactions when the transition probability is small. We use extensively the high 
sensitivity of the NOMAD experiment to illustrate potentially large effects from small, direct flavor 
violation. In the purely leptonic sector, we find that the clean and v e beams from a /i + — 
collider could provide the sharpest tests of direct flavor violation. 

I. INTRODUCTION 



The Superkamiokande collaboration's zenith angle analysis of its data Q, may sway even a skeptic to the view that 
neutrino oscillations have been observed. The precise form of the neutrino mass matrix and the number of neutrino 
^ ' species remain undetermined, however. The confusing and complicated nature of the whole collection of neutrino data 
i in laboratory and astrophysical settings makes for an exciting playing field for workers trying to establish a "standard 
model" of neutrino physics. Vigorous efforts have been made to determine whether one needs only three flavors of 
neutrinos, mixed either two-by-two 10] or with some form of full three- state mixing Q , Q or whether a fourth, sterile 
neutrino is needed in the mixture g7 @ • There is not yet a compelling phenomenological or theoretical case for any 
of these reasonable, and partially successful, approaches. 

Direct violation in neutrino interactions is a feature of the flavor-mixing puzzle that has received relatively little 
attention. Most models that predict the existence of neutrino mass and mixing, and consequently neutrino oscillations, 
also contain direct neutrino flavor- violating interactions. In fact, some interesting models have no neutrino masses 
at tree level, but flavor- violating interactions generate them in loop graphs 0, II, pfj[ , |yj . Moreover, with the 
Oh, help of resonant enhancement, massless neutrinos can mix in a non-trivial way in certain models and produce flavor 
transitions in flight Jl2| . It seems natural to extend the phenomenological framework for oscillations of propagating 
neutrinos to include the effects of direct interactions in a way that allows one to survey all effects at once. A formalism 
to carry out this program is sketched and illustrated in p3| , and we extend the range of applications and uncover 
several new features in the present paper. Specific models are not discussed here, since we emphasize the model 
independent features of the oscillation-plus-direct flavor violation analysis. Bounds on parameter combinations from 
a given model can be directly obtained from the coefficients of our effective four Fermi Lagrangian. 

As in [^3|, our primary concern is with the accelerator experiments [[hi} |H| |^6| . The experimental constraints 
on muonium-antimuonium conversion are now so tight nj\ that a purely direct interaction explanation of the muon 



decay- at- rest (DAR) signal reported in ][14j can be rather convincingly ruled out 18 1. Here we emphasize effects that 
are significant when combined with the oscillation phenomenon. For example, we show in section IV that direct effects 
can kill oscillations in special circumstances. We also find that in high sensitivity experiments, where tight limits 
in regions of sin 2 26- Am 2 space are achieved, small, direct flavor violation can change a boundary by more than an 
order of magnitude. We illustrate this by showing examples of the effects on the large Am 2 , small sin 2 28 boundary 
of the i/fj, v T mixing set by the NOMAD collaboration . In another sensitive comparison, we show the power of 
comparing the "wrong flavor" appearance signals from the clean and v e beams afforded by a /x + — u~ collider. 

There are several studies where direct flavor violation is considered in the solar and atmospheric cases. In |2fj[| the 
combined oscillation and direct effects are applied to an analysis of the resonant conversion of electron neutrinos to 
other species as the explanation of the solar neutrino deficit plj |22j |23| |24| , while in |^5| and p6| , an explanation 
using resonantly enhanced direct lepton number violation of the zenith angle effect reported in |lj was presented. In 
both cases significant effects were reported, though a complete explanation in terms of direct interactions is probably 
not possible in the solar case. The situation is not settled in the atmospheric case. A critique of alternatives to the 



large v^- v e mixing solutions to the atmospheric neutrino data such as those presented in j25[ , |26| and |27| is given 
in [28] ]. A model using decay that answers the objection in [p8| is described in ]29| . 

We encourage the reader interested in the impact of small, direct flavor violation on the analysis of experiments 
with high sensitivity to " wrong flavor" appearance to go straight to Sec. 4. The background is given in the next 
two sections. In the following section, the parameterization and notation are defined, and the notion of a generalized 
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transition probability factor is explained. In Sec. 3 the experimental constraints on flavor violating parameters defined 
in Sec. 2 are summarized. As just remarked, the formalism is applied to a number of examples drawn from current 
or future accelerator experiments in Sec. 4. We keep the constraints on flavor-violating parameters clearly in mind 
in the discussion of these applications. The formalism and results of its applications are summarized and several 
conclusions are drawn in Sees. 5 and 6. The general forms of the probability factors that apply to the case where fj, 
decay provides the source of neutrinos are given in the Appendix. 



II. FORMALISM 



In this section we develop a compact parameterization of direct interaction effects in neutrino flavor-changing 
processes. We represent the low-energy effective interactions involving neutrinos, charged leptons and first generation 
quarks by the four-Fermi semileptonic (S) and leptonic (L) Lagrangians^ 



C s = 2^G F K h A {kY A P h U ]a v a )[dT A {aP L +(3P R )u^ + h. 



and 



C L = 2V2G F F$'(Jlr A P h U ja v a )(hr A P h ,U rnb v b ) 



t 



(1) 



(2) 



where U{ a is the matrix, unitary in the relativistic limit, that relates flavor as shown below in Eq.(7). Repeated 
indices are summed. The coefficients K and F, whose indices are i,j,k, represent the coupling strengths for the 
different lepton flavor combinations, while the indices a, b, ... label the mass eigenvalues. The coefficients a and 
(3 allow for different strengths for L and R couplings to the quark currents. The Lorentz structure of the bilinear 
forms is labeled by A = S,V or T and Pt denotes the left- and right-helicity projections. The expression in Eq.(2) 
is a generalization of the generic muon-decay, four- Fermi interaction to include lepton flavor violations of all types. 
By a Fierz transformation one can show that F^ L and F RR are both identically zero. Restricting application of 
Eq.(2) to muon decay and (unobserved) massless neutrinos, one can show that it is not possible to test lepton number 
conservation from the available observables |3(J . 

We illustrate the notation by applying it to the standard model (SM) effective, low-energy Lagrangian. The leptonic 
neutral-current term, 



£ L SM C = V2G F [k(2s 2 w ^P R + ^(2s 2 w ~ l)P L )U](T7-^P LV3 ), 
can be Fierz transformed into an equivalent charged current form 



£ L SM C = V2G F [[k{2s 2 w - l)^P LVj ](VJ^P L li) - 2[l l 2s 2 w P L ^]{T7]P R h)}. 
Next we add the SM leptonic charged-current effective Lagrangian 

C^f = 2V2G F {k^P L v l ](l7-^P L l 3 ), 
and the semileptonic effective Lagrangian 

C S SM = 2V2G F [u 1 , l P L d](T 3 ^P LVj ) + h.c. 
The F and K coefficients can now be read off from the SM effective low energy Lagrangian: 



(3) 
(4) 
(5) 
(6) 
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(24-l)/2 





1 The CKM factors multiplying Gf play no direct role in our discussion, so they are suppressed in the notation. 
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A. Lepton Flavor-Changing Transitions 



In the usual analysis of lepton flavor oscillations, the neutrinos are treated as massless in the matrix element 
kinematics. Approximating the plane wave phase factors for the propagating neutrino to leading order in the masses, 
one factors out the transition amplitude to write^J 

a a 

However, this factorization is not valid when the neutrino masses are taken into account [ |3l| . Though we work in the 
ultra-relativistic limit, where the neutrino masses can be ignored in the arguments of the matrix elements, we find 



the process dependence as discussed in 1 31 to be a useful setting for the intermediate stages of our development. We 



generalize the SM weak-process initial states to include new physics; for example, 

W»)WP ~Y,\^^V + \£Im\* + ) ^EaK)(Va,» + \C S \TT+), (8) 
a 

with L s defined in Eq.(l). Similarly, we create a weak- process final state for the detector and use it with Eq.(8) to 
construct the transition amplitude. More generally, denoting initial and final states by \I s ' d ) and |,F s ' d ), we define 

M s a = (F s (v a )\£ s +£ L \I S ), (9) 

and 

Mj = {F d \C s +C L \I d (v a )), (10) 

as the source and detector transition matrix elements involving a mass eigenstate of v a . We can write the full transition 
amplitude, including direct new interactions, from creation to detection of the neutrino as 

(vj(t)\"M)Ni = E^>"> e ~""" t/2 >>^ =M d a e— 2 ^ E M s a , (11) 

a 

where a sum over the label, a, of mass eigenstates is implicit here and in what follows, and NI indicates that new 
interactions are included. 

As an example, let us consider a SM process that is a background to electron appearance experiments^, namely 
7r + — > e + v e followed by v e Ni — > e~Nt, where Nij designate inital and final hadronic (nuclear) states. We expand 
the SM to include the possibility of neutrino oscillation for illustration. The transition matrix elements for the source 
and detector processes are 

M s a = 2^2G F K^ 11 (e+^\(h^P L U la ^\0)(0\O v \7T+) (12) 

and 

M d a = 2^2G F K L ll {N f \0^\N t ){e-\h ltl P L U la v a \v e ), (13) 

where Oy designates the quark current operator cIj^Plu appropriate to the SM in Eq.(l). In the following, we will 
continue to use O to designate the hadronic current operator. In the ultra-relativistic limit, the neutrino masses are set 
equal to zero in the spinors and only the leading phase dependence on masses is kept. The transition amplitude squared 
can then be factorized into a SM product of matrix elements squared times the oscillation probability, involving a 
sum over mass eigenstates: 

\ M d, e -imlt/2E Ms j2 = ( 2 ^2G F ) 4 \(e^P L vy(e lx P L v)(0\O v \ud^ (14) 
where we have used the fact that Kyu = 1 in the SM. Schematically we can write the above as 

\ M d e -imlt/2E M sj2 „ $ ea (y eNi _ e-Nf)P e ^ e . (15) 

In Eq.(15) <E> e designates the v e flux, a(v e Ni — > e~Nf) the SM electron neutrino charged current cross section and 
P e - >e = \Ui a e~' lm <' t / 2E Ui a \ 2 designates the probability that an electron neutrino produced at the source appears as an 
electron neutrino at the target. Our next task is to achieve an equally transparent factorization that includes direct 
violation of lepton flavor. 



2 In the ultra-relativistic limit, one may take t=L, the propagation length, in the following expressions. 

3 These are experiments that look for signals of neutrino flavor j interactions in a beam of neutrinos created with flavor i 7^ j. 
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B. Lepton Flavor-Changing Probability Factor 



To identify a flavor-changing factor that includes direct flavor violation effects and plays the same role as the 
oscillation probability factor P e ^ e in Eq.(15), we use the amplitudes defined in Eqs.(9), (10) and (11) and follow the 
pattern of our preceding discussion of Eqs.(12) - (15). Our new "probability" factor should reduce to the oscillation 
probability when the direct flavor-violating couplings are turned off, of course. For this purpose, we continue to use 
the pion decay source and nucleon/nucleus detector example, but now we will extract a lepton appearance probability 
factor from combined oscillations and direct interactions. We adopt an illustrative model with a V-A current with 
Ky 2l =^ 0, with I — t or e but with all other non-SM coefficients K^=0. This model produces a ir + — > /i + + vi 
transition at the source but leaves only the SM process vi + Ni — > l~ + Nf active at the target. We shall refer to this 
as a "source-only" situation. The electron production rate at the detector is proportional to 

\ M d e -imlt/2E M , a \2 = \(l-Nf\£ s \v a N i )e- im « t / 2E {fi+v a \L s \ir+)\ 2 , (16) 

where the right-hand side, after spelling out the terms from Eq.(l) and taking the ultra-relativistic limit, reads 

(2y2G F ) 4 |7 7A F L ^<iV / |^|iV l )| 2 |(0|^|^)l7 > P LAi |2| ([ / 2a + K L 2iUla)e - lm lt/2E Uta?i (17) 

where Kyu = 1 has been used. The v\ appearance probability factor, including the direct flavor- violation coefficient 
Ky 2l , is given by the final expression within absolute magnitude signs in Eq.(17). The first term in this expression is 
the usual oscillation amplitude for transition from to vi between source and detector. 
To make the notation more compact in the following discussion, we define 

(m 2 - ml)t/2E = 2x 

K^ 2l = tamPe 2l,p , 

where the choices of a and b in the definition of x and the values of ip and cj> depend upon the situation, and we will 
illustrate the two-flavor mixing and the three-flavor mixing with dominance by one mass scale || . We have chosen the 
parametrization of Ky 2l by tamp to emphasize that the flavor-violating amplitude could, in principle, be larger than 
the flavor-conserving one. In practice it is constrained by experiment to be small compared to the standard model 
amplitude, but the parametrization reminds us that the direct flavor- violating amplitude cannot be represented as 
a unitary rotation to a "source basis" except in the small tanip limit with (p negligible. In addition, the special 
circumstance must occur that the effective flavor-violating four-Fermi interactions have a V-A structure to match 
the SM structure. Nonetheless, we will refer to the transition factor modulus squared as a probability or probability 
factor in what follows. As we show below, the survey of direct effects on the standard two- and three- family mixing 
plots with Am 2 vs. sin 2 (29) contours or tan 2 9i^ vs. tan 2 &23 f° r fixed 5m 2 , for example, can be straightforwardly 
extended to include direct lepton number violation. 
Isolating the probability factor in Eq.(17), we write 

P^i = \U 2a e- lm « t/2E U; a + tan^e 2 ^U lb e~ m ^ 2E U; b \ 2 , (18) 

no sum on I, which labels e or r. With <fi = 0, in "mock" unitary formQ we have 

P^ t = \cosi>U 2a e- m ^l 2E Ut a + S inmbe- v<t/2E Ut b \ 2 /{co S ^) 2 , (19) 

and, approximating the overall (cosip)~ 2 factor by 1 to order ip 2 , Eq.(19) can also be recast as 

P^i * \V 2 2U 2a e- im ^ 2E U^ a + V 2l U lb e- im W 2E Uf b \ 2 . (20) 

We have identified V22 = cosip and V2i=sinip in Eq.(20), which suggests the interpretation that the /i from n decay 
is accompanied by a with amplitude V22 and followed by propagation and an oscillation to v e with amplitude 
U2 a e~ m2at ^ 2E Ui a or is accompanied by a v ; / 7^ fx with amplitude V21 and then followed by v\ propagation with ampli- 
tude Ui b e~ lmbt ' 2E U* b to remain vi. This picture, though only approximately valid and only then in the circumstances 
described above, is useful for seeing how the direct lepton flavor violations can be worked in with the oscillation 
between flavors in a reasonably seamless fashion. 



4 For present purposes, we assume <j> is constrained to be small by experimental limits on CP violation. 
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III. CONSTRAINTS ON THE STRENGTHS OF DIRECT LEPTON-FLAVOR VIOLATIONS 



The identities of electron neutrino, muon neutrino and tau neutrino are preserved to a good accuracy in decays and 
collisions. Within neutrino data alone, the constraints on the sizes of lepton-flavor violations, though tight, often leave 
room for violations at the 1/2% level in cross-sections and rates. In the realm of purely charged lepton processes, 
on the other hand, the high degree of experimental control has lead to fantastically small limits on some of the 
ratios of lepton-number- violating to lepton-number-preserving interaction strengths. In this section we briefly survey 
the constraints that affect our analysis the most directly. The considerations presented here govern our choice of 
parameter values in the following sections as wc illustrate some of the possible effects that direct interactions produce 
in concert with oscillations. 

Weak SU(2) invariance of any new physics interaction Lagrangian that produces lepton number violation will 
generally relate the purely leptonic processes to those involving neutrinos . The severe experimental limits on the 
kinematically allowed leptonic processes then translate into limits on processes involving leptons and neutrinos, up to 
group theory factors and SU(2)-violating mass-splittings among members of boson multiplets that mediate the lepton 
number violations. The purely leptonic processes whose experimental limits impose the strictest bounds on lepton 
number violation in muon-source experiments are muonium to antimuonium conversion JlTy , /i — * eee, r — > fiee and 
t —>■ /i/ie (3^1 . There are similar constraints (3^] on the strengths of semileptonic processes violating lepton number 
that follow from r — * 7r + e, r — ► 7r + fi and from fi conversion to e on Ti nuclei [ p3| . The constraints are summarized 
in Table 1. The "experimental constraints" shown in Table 1 are actually those that apply to the charged-lepton 
processes listed. The relationship to the bounds on the lepton-number-violating coefficients for the neutrino processes 
is somewhat indirect, since group theory factors, ratios of masses of virtual bosons mediating the processes and effects 
due to differences between, say, S and V structure must be included. Allowing a possible factor of two from Clebsch- 
Gordon coefficients and a generous factor of two in the ratio of masses of exchange bosons within the same SU(2) 
multiplet, the coefficients shown in Table 1 are bounded by roughly eight times the purely charged- lepton process 
limits listed there Jl8[ . This value is given in the column labeled "model independent constraint" . The precise value 
of the factor for a given process is model dependent. In Table 1 and in the rest of the paper, a superscript L is to be 
understood if none is shown explicitely. 



coefficient 


process 


experimental constraint 


model independent constraint 




-FV2111 


/1 — ► eee 


1.0 x 10 _b 


8.0 x 10~ b 




-FV2112 


muonium-antimuonium 


3.0 x 10" 3 


2.4 x 10" 2 




-FV2113 


r — > jj, + ee 


2.9 x 10 -3 


2.3 x 10~ 2 




-FV2213 


r — > e~ fi^i 


3.3 x 10" 3 


2.6 x 10" 2 




-FV2311 


t —> n~ee 


3.2 x 10" 3 


2.5 x 10" 2 




-FV2312 


t — ► e + /ifi 


2.9 x 10" 3 


2.3 x 10~ 2 




Kv21 


H — e conversion 


1.8 x 10" 7 


1.5 x 10~ 6 




Kvn 


r — » e7r° 


8.2 x 10" 3 


6.6 x 10~ 2 




Kv32 


r — > jj,Tv a 


8.5 x 10 -3 


6.8 x 10~ 2 





TABLE I. Limits on charged lepton processes and the corresponding neutrino process coefficients F and K. The model 
independent constraints on the F and K coefficients are taken to be roughly eight times the charged lepton process limits. For 
K constraints we take h=L, a=l and /3=0 in Eq.(l). 
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IV. IMPACT ON ANALYSIS OF APPEARANCE EXPERIMENTS 



The very term "neutrino oscillation" implies path-lcngth-dcpcndcnt variation of the probability that a given flavor 
of neutrino appears in the beam. Moreover, neutrino oscillation and neutrino mass are so tightly linked that evidence 
for the former is considered tantamount to proof of the latter - certainly in the case of vacuum oscillations. Conversely, 
the absence of oscillations in a neutrino flavor-violating effect is tantamount to elimination of neutrino mass as an 
explanation of its origin. This is not necessarily so when flavor violation is expanded to include direct interactions. 
This point is among a number that we make in the present section. The examples chosen are all consistent with the 
bounds described in the previous section and summarized in Table 1. 

At the beginning of this section we take the neutrino source to be ir- decay. We explore the interplay in the 
Vfx <-> v T case between the mass-induced, oscillating amplitude and the directly-induced, non-oscillating amplitude. 
While it is true that small masses generally lead to vacuum oscillations, it is not strictly true that the absence of 
oscillations proves that the neutrinos are massless. In a sense this is a complement to the well-known result that 
massless neutrinos can oscillate as they pass through matter. These effects are shown in Figs. 1. and 2. In Fig. 3., 
we show the result of including the flavor- violation parameter "axis", ip, in the analysis of the probability bounds in 
the sensitive NOMAD experiment. The result is dramatic. Then we make the point that the bound on the oscillation 
mixing angle for fixed Am 2 depends in general on the value of the direct flavor violation parameter ip. This is shown in 
Fig. 4, which includes the possibility that the same direct flavor violation occurs at both the source and the detector. 
This discussion is followed by the expansion of the analysis to the three flavor situation. Figures 5. and 6. show how 
different tan 2 9is - tan 2 #23 boundaries appear as different fixed ij) and Am 2 "slices" of the parameter space are taken. 
Figure 5. is appropriate to the upcoming MiniBooNe experiment, while Fig. 6. applies to the reported NOMAD 
probability bound. 

In the last part of this section, we show the power of the clean and v e beams from proposed storage rings 

to make sensitive tests for direct flavor violation in the purely leptonic sector. Of particular note is the prospect of 
advancing another order of magnitude into the parameter space of direct v e <-» flavor violations. 



A. 7r-decay as the Neutrino Source 



Restricting ourselves at first to mixing of two mass eigenstates, we write Eq.(19) for the probability of lepton 

V21 



appearance in the pion decay as the source in an "all angles" form and defining Fy 2l = tanipe 21. 



P^i = \cos^{-cos9sin9e lx + sin9cos9e~ rx ) + sin^e 2li '{cos 2 6e ix + sin 2 9e- ix )\ 2 co S - 2 %p. (21) 
As defined above, x = (m 2 — m\)t/AE in applicaton to Eq.(21), and the two-flavor mixing matrix is written as 

U=( C0S ° a Sin 2V (22) 
V — sin 9 cos 9 J v ' 

Rearranging terms and consolidating them, we arrive at the rather transparent form of Eq.(21) 

2 sin29sin2(9 — ip)sin 2 x 
Pa^i = tan ip H — h 4tanipsin29sinq>sinx(cos29sin(psinx — cos<pcosx). (23) 

COS Z %! 

Equation (23) has the obvious and expected feature that if ip=Q, we have Pu-^i = sin 2 29sin 2 x, the usual two-flavor, 
pure mass- mixing, oscillation formula in terms of the mixing angle 9 and the factor x = Am 2 L / 4E . Equally obvious 
and expected is the relationship P^i = tan 2 ip that holds when 9 = or x = 0. What is not expected is that, when 
sincp = 0, P^i = tan 2 %p, independent of x, when 9 — ip = irn/2, where n is an integer. Looking back at Eq.(21), we see 
that when 9 — tp — irn/2 and when </>=0 the coefficient of e %x from the — » vi oscillation cancels against its coefficient 
from the v\ — > v\ term in the direct flavor violation amplitude. The remaining overall phase from the e~ %x factor 
disappears in the modulus squared and one is left with simply P^i = tan 2 ip = tan 2 9. This somewhat surprising 
result in the case where there is direct flavor violation in the DIF source dramatizes the implications of Eq.(23) for 
interpreting signals for oscillation, or lack thereof, in variable baseline experiments. There are counterparts to this 
source effect in the muon decay case JlJ] |l{J as well as in the cases where direct flavor violation occurs only at the 
detector or in both the detector and the source. We will comment further on these situations and on three-flavor 
mixing below. 

The condition for exact cancellation of the L/E dependence is unlikely, of course, but the interplay between the tp 
and 9 dependence is generic, and it affects, possibly radically if mixing angles are small, the interpretation of signals 



G 



that show variation with path length. We illustrate the L/E-dependent effects caused by the interference between the 
pure oscillation term and the direct flavor violation in Figs. 1 and 2. 



0.007 



0.006 



psi=0, phi=0 
psi=-0.02, phi=0 
psi=-0.02, phi=0.1 
psi=-0.02, phi=-0.1 
psi=-0.01, phi=0 
psi=+0.01, phi=0 
psi=+0.02, phi=0 




FIG. 1. Dependence of P M _+; on x for 8 = 0.03 and several different values of ip and with Am 2 = lev 2 . A Gaussian smearing 
model has been adopted with a = 0.016, appropriate to the NOMAD experiment, which we will use as an illustrative example. 
The top three curves, which are barely distinguishable from one another, show the very weak dependence on the CP-violating 
phase cj> 



Figure 1 illustrates the variation of P^i with L/E for fixed 9 = 0.03, or sin 2 29 = 0.0036, for various values of ip. 
Figure 2 shows the variation for ip = —0.02 and various values of 9. Several combinations of 6 and ip can lead to a 
given curve, which suggests that it would not be straightforward to disentangle the oscillation parameters from the 
comparison of the x-dependence of the probability with data if small direct effects were included in the analysis. 



Conversion Probability vs b with psi=-0.02 
0.045 | 1 1 1 1 




0.005 



1 ^ ' ' 1 J '- - ' 

0.5 1 1.5 2 2.5 3 

b=UE 

FIG. 2. Dependence of P M ^; on x for ip — —0.02 and several different values of 6 and with Am 2 = lev 2 . 



For completeness, we also show in Fig. 1 the (weak) influence of 4> on the two-flavor oscillation amplitude. Com- 
paring the top three curves, we see that including the ^-dependent effects for reasonable <j) values leaves the picture 
essentially unchanged, as one expects when a small CP-violating phase rests upon a small flavor-violating amplitude. 
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A CP-violating phase of order one is needed if it is to make an observable impact. We reserve the analysis of effects 
from large CP-violating phases |3 for a future study. 

We see that the behavior of the probability is governed essentially by the second term in Eq. (23). This term makes 
it clear that even small ip values can have a large effect if the mixing angle 9 is of the same order of magnitude as ip. 

Are there circumstances where the sizes of these direct effects could be big enough to be observable? A glance at 
Table 1 shows that Ky2i (=tanip) is constrained to be far too small, even with generous allowance for SU(2) breaking 
effects, to modify the pure oscillation picture of from ~k decay oscillating to v e . The Kv32 bound is much looser, 
however, and interesting effects could occur in — > v T two-flavor mixing with tamp < 0.02 if # MT were of the order of 
0.1 or less. A large mixing angle is required between and another species by the two-flavor fit to the atmospheric 
neutrino anomaly. Thus mixing with a sterile neutrino is required, in order to have a direct interaction effect that 
shows up in a (weakly) mixed za, <-> v T sector and significantly modifies the v T appearance signal from a pion decay 
source of v^. For small appearance probabilities, the change in signal depends rather sensitively on the value of the 
direct interaction strength. We give some detail in an example in the /i decay case in subsection B below. 



1. Impact on Boundaries in /i — r Mixing Space 



Failure to detect v T 's in an appearance search allows one to set confidence level curves in /i-t mixing space, and this 
gives another slant on the application of our formalism. The recent NOMAD results give the smallest probability, 
and tightest large Am 2 limits on the allowed region of parameter space for two-flavor, <-> v T mixing. To illustrate 
the impact of small direct interaction effects on the NOMAD bound, we approximate their boundary curve, which 
corresponds to a P=0.0006 appearance probability, by the simple Gaussian smearing model |32| with parameters fit 
to reproduce the main features of the NOMAD boundary in sin 2 8 vs. Am 2 plane for mixing of and v T . Keeping 
these parameters fixed, we replot the contour in the sin 2 28-Am 2 plane for several small values of ip. As before the 
small <j) effects are not interesting for our present point, and we set 0=0. The result is shown in Fig. 3. 

NOMAD Mu Neutrino->Tau Neutrino Limits, P=0.0006 
3 I : — i 1 1 q n 1 — i 1 1 1 

: i psi=0 

I ! psi=0.02 

I psi=0.01 

2.5 - I psi=-0.02 - 



1.5 



O 
O 



0.5 




- 



. .5 I 1 1 1 1 1 1 1 1 1 1 

-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 
LOG(sin 2theta-squared) 

FIG. 3. Model of the NOMAD 90% C.L. limit on — * v T oscillations for several values of ip with <f>=0. The curves correspond 
to P=0.0006. 

The boundaries of the null search results are significantly revised even when small direct effects, consistent with 
the bounds from other experiments, are introduced. NOMAD's high sensitivity to small appearance probability in 
the large Am 2 region is the reason that the inclusion of small direct amplitudes has such a pronounced effect. 

A few comments on the shape of the different curves is in order. NOMAD has a small average L/E value, so 
it is sensitive only to large Am 2 . Their L/E distribution is broad, so the smearing almost completely damps the 
oscillations in the boundary curve at large Am 2 . The Gaussian model is even more extreme in this respect than the 
actual fit to the data. Because of the factor sin2(8 — ip) in the expression for P^— > T > for < 9 < tt/2 the positive ip 
values reduce the amplitude and the sin 2 x factor has to compensate. This eventually drives the curve to higher Ato 2 
for a given 6 value. The opposite behavior occurs when ip is negative. The role of the term tan 2 ip is most apparant at 
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large values of 9, where the ip in the argument of sine is inconsequential, but the addition of tan 2 ip to the probability 
allows a fit with somewhat smaller Am 2 at 9 = ir/2. The crossing of positive ip curves with the ip =0 curve is forced 
by the large Am 2 and small Am 2 behaviors just described. 



2. Direct Flavor-Change Bounds in ip — 6 Space 



Proposed high-sensitivity experiments to probe smaller Am 2 and sin 2 29 regions can also place discovery limits and 
upper bounds on direct, flavor- violating interactions involving neutrinos. For example, Fermilab J35|] and CERN |3^] 
proposals aim to push P^ e bounds down to 10~ 4 at 90% C.L., while the Fermilab-Soudan experiment, MINOS [p8[ , 
is shooting for a bound of 10~ 2 on P^ T . Similar sensitivity is proposed in muon collider sources of pure v e and 
beams. 

Let us consider the situation where z/^'s originate from 7r decays, which will be the case in the MiniBooNE experiment 
p5j . If the new flavor physics is only at the source, Eq.(23) applies. If one has the same new physics amplitude at 
the source and detector, the corresponding expression is 



P„_>; = Atan 2 ipcos 2 x + sin 2 28sin 2 x — 2tamfjsin29sin2(f>sin2x. 



(24) 



Assuming that the CP-violating phase, <fi, is small, we can readily illustrate the influence of the x and 9 values on 
the limits on ip imposed by a given 90% C.L. bound on P^_>;. In Fig. 4 we show the boundaries in if) — space 
for a fixed value of x=1.27(L/E)Am 2 =l when P M _»i < 10~ 4 is imposed. Fig 4 gives the source-only case and the 
source-plus-detector case boundaries. The specifications of MiniBooNE are L ~500m and typically 0.5<E<1.0 GeV 
p5[, so x~ Am 2 for purposes of translating the graphs to MiniBooNE's capabilities. 
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FIG. 4. Contours of 90%C.L. upper bounds imp — 6 space for x=l and P M -» e =10 4 . The cases where the new physics is at 
the source only and where it is at both the source and detector with equal amplitudes are shown. 



In Fig. 4 the areas within the "egg-shaped" regions are allowed at the 90%C.L., and the maximum allowed direct 
flavor violation parameter is evidently a function of both 9 and x. It is generally assumed that the best bound on 
the flavor-violating amplitudes is obtained when the oscillations are not present, but the "tilted egg" in Fig. 4 shows 
that this is not necessarily the case. The interference between the oscillation and direct amplitudes in the source- 
only situation makes the biggest direct flavor violation effect occur at non-zero 9. When the same flavor violation 
is assumed to apply at the source and the detector, the symmetry of the set-up ensures that the maximum allowed 
value of ip occurs at 9=0, which is the usual expectation. 

With the Am 2 value chosen, the bound \ip\ <1.6xl0 -2 results, which is an order of magnitude stronger than the 
current direct bound from neutrino processes, but still much weaker than the bound that can be inferred from the 
limits on (i — > e transitions in heavy nuclei as listed in Table I. 
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3. The Three-Flavor System 



Distinct new possibilities are created when e, [i and r-flavors are all included in the picture with mixing plus direct 
flavor violation. Here we re-examine the 7r-source, e and r appearance possibilities with three flavors, using the "one 
mass-scale-dominance" model |3). 

Since the experimental constraint on Kyi\ makes it irrelevant for our purposes (see Table 1), we do not include 
it in the expressions below. With flavor violation at the source, and restricting ourselves to V-A structure, we need 
only consider the Ky%i and Kv23 coefficients. To the first order in flavor violation, we have 

P^ e = sin 2 28i3sin 2 (923 + ip^sir^xcos -2 ^^ (25) 

and 

Pfi-^r = tan 2 ^^ + Asin 2 xcos 2 9i3Cos623Cos~ 2 ip fJi sin(923 + ipfj,)(—sinip^ + cos 2 9 \3Cos9 23 sin{9 23 + Vv))> (26) 

where K23 = tanipfi and the mixing-angle convention is that of p2| ; namely U\3=sin9i3e~ Sl3 , [^23= sin923Cos9\3 and 
U33—cos9i3Cos923- We have assumed <5i3=0 and real Kytj values. P^ e (P^ T ) includes the amplitude that is 
produced at the source and oscillates to v e {v r ) plus the amplitude that v T is produced at the source and oscillates to 
v e (remains v r ). Note that, as in the two-flavor case, the special condition 9+ip=mt kills the oscillation amplitudes, 
and in the /x — > r case there is a second condition when this can happen. The lack of symmetry between P^—, e and 
P^r results from our neglect of the K21 coefficient in the amplitudes. Figure 5 shows the effect of choosing different 
x values (i.e. different Am 2 values for fixed L/E) and non-zero tp^ values on the P Ai ^ e =10 -4 boundary in the tan 2 9\3 
vs.taro 2 (?23 plane, which is appropriate for the MiniBooNE [B5j parameters. 




-5 1 1 1 1 1 1 1 

-6 -4 -2 LOG tan«S(theta23) 2 4 6 

FIG. 5. Contours of 90%C.L. upper bounds in tan 2 8is vs. tan 2 023 space for several sets of tp^, Am 2 values and with fixed 
value P^ e = 10- 4 . The new physics effects are taken to be at the source only. 

For Am 2 = ley 2 the change in the small tan 2 #23 region as Vv is turned on is noticeable already at ip mu =0.01. The 
distortion grows rapidly with ip, and the LOGtan 2 923 minimum value moves well below 10 -6 when -0^=0.02 and 
Am 2 =l. The symmetry of the graph about the tan 2 9i3=l line is a consequence of the equivalence of the 9 values 
below and above 7r/4 caused by the sin 2 29 factor. 

The influence of flavor violations in r-flavor processes is not so tightly constrained experimentally as for strictly 
e-plus-/x, and we look at this situation next. We contrast the bounds on r-appearance probability in the cases that 
the flavor violation is at the source alone and the source and detector both. In Fig. 6 we again look at the tav?9\3 
and tan 2 923 plane boundaries for P M ^ T =6 x 10 -4 , as in the NOMAD experimental boundary discussed in IV. A. 1 
above in the two-flavor model. Plots with Ato 2 =9,6 and 3eV 2 are shown, with oscillation only in the 9 and 6 eV 2 
cases but with ^=0.02,0.0 and -0.02 for the source-only case at Am 2 =3eF 2 and, for comparison, ip^ =0.01 and -0.01 
in the source-plus-detector case at Am 2 =3eV 2 . 
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FIG. 6. Contours of 90%C.L. upper bounds in tan 9m vs. tan 823 space for Am 2 =9,6 and 3eV 2 . The curves are designated 
as source or source and detector according to whether flavor violation is included at the source only or both the source and the 
detector. In parethesis the first number is the Am 2 in eV 2 and the second is the value of Vv- F° r ^=0 there is no direct flavor 
violation, so the "source" and "detector" labels are not meaningful. P M ^ T =0.0006, appropriate for NOMAD. 

The most noteworthy feature of Fig. 6 is that there is a "confusion effect". The boundary curves for Am 2 =3 and 6 
eV 2 for different flavor violation situations crowd together and give a portrayal of the complications that arise when 
direct flavor violations and oscillations become competitive. The smaller values of ip used in the "source-plus-detector" 
situation compared to the "source-only" are forced by the lack of solutions to the probability equation when \tp\ is 
0.02. The reason is that there is a leading term 4 x tan 2 ip, so that there is no way to obtain P=0.0006 when the lead 
term is 0.0016 and the sin 2 x term is as small as it is in the NOMAD experiment because of the small L/E. 



B. /x-Decay as the Neutrino Source 



With muon decay as the neutrino source, there are electron and muon neutrino appearance examples within this 
two-state system involving the flavor-violating coefficients -FV2111, -FV2212 and i*V2ii2- We give general expressions for 
cases where u e — ► v^, v e — > v T , v ^ — > e , and — > v T in Appendix A. Again restricting ourselves to a V-A form for the 
new interactions, the general form of the probability factor for a — > e transition is 

P^ e = I — 2isinxe~ ix sin9cos9Fv22ii + (1 - 2isinxe~ lx sin 2 6) Fv2iu\ 2 + 

I - 2isinxe~ lx sin9cos8 F2212 + (1 - 2« sinxe~ lx sin 2 9)F v 2112I 2 • (27) 

The corresponding expression for P e ^n is obtained by the interchanges -FV2111 -FV2212 and cosO <-» sin6. As in 
the 7T-decay case, Fig.(l), unless phases of the F's are of order 1 they play an insignificant role in the probability 
factors and we take them to be zero in the following discussion. Referring to Table 1, we see that we may drop the 
Fv2in=tcLn'ip ee term since it is several orders of magnitude smaller than the other lepton number violating coefficients 
in /i-decay. We take Fv22ii=l an d expand the resulting expressions using the parameterizations Fv2U2=tanip elJ , and 
Fv22i2=tanip flfl . The resulting expressions are given by 

P^ e = tan 2 %l) efl + sin 2 xsin 2 29(l + tan 2 ^^ — tan 2 ip etl — 2coi26*tara/> wt tara/> e /j) (28) 

and 

P e ^fi = tan ip^ + (sin 2 xsin29sin2(9 + ip^^cos^ 2 ^^^ + tan 2 tp etl (l — sin 2 x sin 2 29)). (29) 

Dropping the ipee in the expressions for P^ e and P e ^^ produces the lack of symmetry between the two equations. 

It is clear from Eq.(27) and the restriction I-FV2112I = \tanip ei i\ < 8x0.003 =0.024 that the effect of direct interactions 
on the amplitude of a — > e oscillation will be important only if it turns out that P^ e < 5xl0~ 4 and experiments 
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can explore that region |37J. We return to this issue below. The e — > fi oscillation amplitude can be more strongly 
affected, because the branching fraction for fj, — > e + + T7^ is not directly constrained by bounds on isospin related, 
purely charged-lepton processes. Therefore, on these grounds alone, sizeable effects cannot be excluded in the v e — > 
appearance case.0 The L/E dependence of P e ^^ in Fig. 7 demonstrates the strong effect that the ip value has on the 
amplitude. 



0.08 



Probability of e to mu transition, mu decay 




FIG. 7. Dependence of P e ^ M for a /x-decay source on L/E for 9 = 0.1, ij} eil = 0.0 and several different values of t/j^ and with 
Am 2 = leV 2 . 

Here we show the oscillation probability as a function of L/E for several choices of ip^ with ip e/Jl =0 and 6=0.1. The 
effects of ip ejl are negligible with the chosen values of and 9 in this case also, as one anticipates from inspection of 
Eq.(29) and the bound on tamper mentioned above. Note that the factor sin2{6 + tp ^) in Eq.(29) kills the oscillations 
in the case ^=-0.1, simply because the choice 0=0.01 is made for the graph. 



1. Comparison of is e Appearance to Appearance 



The behavior of P e ^ M shown in Fig. 7 can translate into significant differences in the appearance probability for 
e compared to /i in neutrino experiments whose beams are extracted from muon collider storage rings, for example. 
In Fig. 8 we show the ratio of P^ e . to P e ^^ as a function of tp^, with 0=0.003 and Am 2 =leV 2 , chosen to be in a 
range allowed by LSND and not excluded by other experiments. 



5 In specific models, indirect constraints on -FV2212 
includes a tamp^ ■ 



tampan from jj, — > ey limits apply to a combination of amplitudes that 
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Ratio of Electron Appearance to Muon Appearance 
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FIG. 8. Dependence of the ratio P M -> e / Pe->n for a \i- decay source on ij)^ for sm 2 2#=0.003, and with Am 2 = lev 2 . See 
text for treatment of L/E. 



The parameters chosen for Fig. 8 are guided by a recent proposal for a medium baseline appearance search with v 
beams from [i decay pffl. Below we summarize the relevant parameter values. 



siri 2 26=0 


(L/E) 


a 


sin z 29=l 


events=11600 fi~ 
T>1 events= 5070 e + 


0.66 
0.75 


0.20 
0.26 


v e events=6400 e~ 
VJl events=3280 ji + 



As Fig. 8 shows, the value of the ratio changes rapidly as a function of ip, and offers a possible method to directly 
constrain ip down to 0.01 or less by comparing the — > v e oscillation signal to its inverse. We can put this another 
way by comparing estimates of the number of e~ events with the number of fi + events detected downstream from 
fi~'s decaying in flight. In a purely oscillation picture, the e~ result from oscillating to v e , while the /i + result 
from 171 oscillating to TTJl. In Fig. 9 we show a plot of the number of events expected vs. ip for the sin 2 28 and Am 2 
values assumed, given the SM event rates estimated in the search experiment proposed in |3^| . 



Number of Electron Events and Positive Muon Events 



50 - 
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0.05 0.1 0.15 0.2 
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FIG. 9. The numbers of e~ and fi + appearance events as a function of V'mmi ^^=0, for sm 2 2#=0.003, and with Am 2 = leV 2 . 
See text for treatment of L/E. 

Figure 9. makes evident that interference between the oscillation and direct amplitudes greatly enhances the /i + 
appearance rate estimate as ip grows, while there is little effect (a factor sec^ip^) on the e~ rate. The numbers 
for -0=0 correspond to the pure oscillation numbers given in p7| , but decreased by the factor sin 2 29= 0.003 for 
Am 2 = leu 2 , since the numbers there refer to sin 2 26>=l. We see again that looking for new physics may be quite 
fruitful in the comparisons of appearance signals in "pure" neutrino beams provided by muon colliders. The current 
direct limit on ip^ can be extracted from the limit ]39| <r(7^7 + e~ — ► /i~ +77)/<j(i^ + e~ — * n^+v) = tan 2 ip< 0.05, or 
tanip < 0.22. As Figs. 8 and 9 indicate, finding N e - > N fl + is sufficient to improve the bound roughly to ip < 0.02, an 
order of magnitude better than the bound inferred from Q. Allowing \ip e fi\= 0.025 has little effect on this statement. 
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2. Impact on Boundaries in e-fi Mixing Space 



Referring next to Eq.(28), let us consider what a proposed reach to < 3.5 x 10~ 4 at 90% C.L. J3?J affects 

the impact of direct interactions on the Am 2 vs. sin 2 29 contour. Specifically, how do fi~~ — * e^v^V^ and — > e~v e 17^ 
sources affect v e appearance from /i~ sources (or Z^T appearance from fi + sources). From (28) one can show that 
tan 2 ipefj. 5: Pfi— >e/(l — tan 2 il) ^^sin 2 x) , which essentially means that tan 2 ip efl < P^ e for small tanip^^. In particular 
for Itamp^^l < 0.2 and P^ e =3.5 x 10 -4 , 0.0187 < (tanip efi ) max < 0.0193. This range is slightly less than the model- 
independent estimate, albeit on the generous side, of the bound \ip e ^\ < 0.024 that follows from the experimental limit 
on muonium-antimuonium transition |l7f . 



We show the Pu_> e = 3.5 x 10 4 contour in Fig. 10 with several sets of values of i/j^ and ip, 



a 
o 




-3 5 -3 -2 5 

LOG(Sin A 2 2theta) 

FIG. 10. P M ^ e =3.5 x 10~ 4 boundaries at 90% C.L. for several sets of Vv^i "^m values as shown on the legend. 

The lack of sensitivity to i/j^ with ip ef ^ ~ 0, already evident from Eq.(28) and Fig. 9 is indicated in Fig. 10 by 
the pile-up of curves with Tp efJl < 0.01. The situation changes drastically when ip — 0.018, near its upper bound for 
solutions for P = 3.5 x 10~ 4 to exist for some Am 2 , range of values. As noted in the introduction, possible direct 
interaction contributions are too small to affect the LSND oscillation signal fits. If LSND results are not confirmed 



and the limits drop to the level of Pn 



10 , then the interpretation of those bounds should include the possible 



range of direct interaction strengths allowed by limits from other experiments. Correlated limits on Am 2 , sin 2 20 and 
tan8w can then be studied with the new data. 



V. SUMMARY OF RESULTS 



Direct lepton number violating interactions and neutrino mass terms often go hand-in-hand in physics beyond the 
standard model. With this in mind we developed a general parametrization of lepton number violating interactions 
that fits smoothly with the usual description of neutrino oscillations in terms of mixing angles. We defined generalized 
"probability factors" and illustrated with a number of examples drawn from accelerator appearance experimental set- 
ups. We concentrated on the case where the structure of the effective four-fermion, charged current, flavor-violating 
operators is V-A. This is the cleanest situation and lends itself to a parameterization in terms of angles that makes 
the role of new interactions in the "wrong-flavor" appearance probability factors rather transparent. For example, in 
Eq.(23) the direct effects are expressed in terms of a leading term that gives all of the appearance probability when 
there is no oscillation, a second term that gives the interplay between the oscillation and direct effects and shows 
that there is a special condition d-ip =integerX7r/2 where the usual oscillation term can be zero even when there is 
a neutrino mass difference and a mixing between flavors. The role of the CP-violating phase in the flavor- violating 
coupling is isolated in the last term, which is zero when the phase is a multiple of tt, so the entire expression can be 
decomposed into readily interpretable pieces. 
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Because the limits on lepton flavor violation are generally tight, the impact on the oscillation picture of appearance 
and disappearance is limited to cases where the oscillation probabilities are small and sensitive to the precise values of 
the mixing angles. We illustrate such situations in detail in Sec. 4, where the sensitive dependence of the appearance 
probability on the relative values of the flavor- violating angle ip and the mass- mixing angle are shown in Figs. 1, 2 
and 7 for the amplitude of the oscillation behavior, are shown in Fig. 3 for the the NOMAD, [i—*t, two-flavor mixing 
boundaries and are shown in Fig. 6 for the corresponding three-flavor mixing boundaries in the one-mass-difference 
dominance model. The NOMAD experiment examples are carefully done within the limits imposed on the size of 
the flavor violations, as summarized in Table 1. We present similar considerations for the case where the source is 
/z-decay, and the boundaries appropriate to a proposed y, storage ring environment are shown in Fig. 10, again for 
realistic constraints on the flavor-violation strength. 

Fig. 4 illustrates another new effect that shows that care must be taken in interpreting limits set by wrong-flavor 
appearance searches. The behavior of the P M ^ e = lCP 4 boundaries for source-only and source-plus-detector are quite 
different as viewed in the Q-\j) plane. The former boundary is "tilted" and so the bound on set by the 90% C.L. 
boundary is correlated with the value of so that, contrary to the usual expectation, the bound on i/j is not obtained 
by setting 9=0. The boundary for the latter situation, however, does conform to the usual expectation. 

We show how the comparison of electron appearance and muon appearance experiments in a clean /z storage ring 
environment could be used to increase the sensitivity to direct, neutrino flavor violations by an order of magnitude in 
the discussion of Figs 8 and 9. These show the rather strong sensitivity of the muon appearance effects on the value 
of the violation parameter ip^. 



A theme that recurs throught the analysis in Sec. 4 is that small, direct flavor violation can seriously complicate the 
picture of wrong-flavor appearance in a number of experimentally realistic situations. In short, the interpretation of 
a signal can be quite ambiguous. To sort out the complete picture, a number of measurements at various values of 
L/E in a variety of different channels is needed. Perhaps the cleanest and most flexible environment for such studies 
is provided by a n + (i~ collider with an associated facility for neutrino beams . Studies such as those exemplified 
in | p7[ can put much more stringent tests on direct neutrino flavor violation than currently exist, while at the same 
time allowing detailed oscillation analysis. Direct flavor-violation effects in semileptonic processes require pion and 
kaon beams of course, and the BooNE experiment will provide exploration of a wide range of parameter space, for 
example. 

Though the direct flavor-violating strengths are constrained to be small, the matter-enhancement effects can lead 
to large transition probabilities, as analyzed in |^5) and pq , and for small probabilities the "vacuum" effects in 
accelerator experiments can lead to important modification of the oscillation description. We conclude that the 
crucial role neutrinos play in our understanding of particle and astrophysics requires that data be analyzed with the 
relevant flavor violation included to properly interpret current and future experiments. The present work provides a 
framework for that task and provides vivid, realistic examples of its application. 
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We gather the general formulas for appearance probabilities that apply when muon decay is the source of neutrinos. 
The presence of two neutrinos in the final state, only one of which is detected by a choice of lepton-flavor sensitive 
detector, makes the analysis slightly different from the situation when meson decay provides the source of neutrinos. 



VI. CONCLUDING REMARKS 



APPENDIX A: MUON DECAY APPEARANCE PROBABILITIES 
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P^ T = J2 \Fi jlk U jc e- iE ^U 3c \ 2 . (A4) 
k 

In all of the expressions, the repeated indices within the absolute-square are understood to be summed. By choosing 
a given process, keeping the dominant flavor-transition terms and parameterizing the F amplitudes by angles, the fi 
decay source expressions in the text can be reproduced from these expressions. 
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